Optimal Control at the Quantum Speed Limit 
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Optimal control theory is a promising candidate for a drastic improvement of the performance of 
quantum information tasks. We explore its ultimate limit in paradigmatic cases, and demonstrate 
that it coincides with the maximum speed limit allowed by quantum evolution. 

PACS numbers: 



Engineering a suitable Hamiltonian that evolves a 
given quantum system into a selected target state has ac- 
quired special relevance after the recent advent of quan- 
tum information science [H. Here, the challenge is to per- 
form quantum tasks (e.g. apply a quantum gate) in an 
accurate way while fulfilling the stringent requirements of 
fault tolerance. In this context quantum Optimal Control 
(OC) is considered a very promising tool, and different 
algorithms have been designed with this aim [3, 0] . One 
of them exploits the Krotov algorithm 0, a numerical 
recursive method which seeks the OC pulses necessary 
to implement the required transformation by solving a 
Lagrange multiplier problem [J. This technique has al- 
ready been applied with success to a wide range of quan- 
tum systems [j, 0] . One issue that is not yet fully under- 
stood however is what its limits are, and how these limits 
may be approached. In this work we will show that the 
effectiveness of the Krotov algorithm for quantum OC 
is related to fundamental bounds that affect the maxi- 
mum speed at which a quantum system can evolve in its 
Hilbert space. Besides being of interest from a theoreti- 
cal perspective, the discovery of such a constraint is also 
important for practical implementations. 

For time-dependent Hamiltonians, bounds that relate 
the transition probabilities of a quantum system to its 
mean energy spread were set by Pfeifer [6|] and Bhat- 
tacharyya [7|, more than fifteen years ago. For time- 
independent Hamiltonians, these results have been ex- 
tended to include dynamical constraints that involve also 
the energy expectation value of the evolving system Q. 
Moreover in a specific case Khaneja et al. evaluated the 
minimum time required to implement a given quantum 
transformation In the light of these results our aim 
is to explore the very limit of OC. In particular, we are 
interested to see whether the Krotov algorithm [2] allows 
one to attain the ultimate bound set by quantum me- 
chanics (for which we borrow from Q the term Quantum 
Speed Limit (QSL)). 

Several attempts to reconcile accuracy and speed in 
quantum control have been proposed so far (see [lol | 



and references therein). In particular, Carlini et al. Q 
cast the time-OC problem into the commonly termed 
quantum brachistochrone problem: exploiting the varia- 
tional principle they produce a collection of coupled non- 
linear equations whose solution (when it exists) yields 
the required optimal time-dependent Hamiltonian that 
minimizes the time evolution while satisfying certain con- 
straints on the available resources. Our approach differs 
from that of Ref. Q since we do not treat the duration 
of the process as a variable that enters in the optimiza- 
tion process. Instead, we set it to some fixed value T and 
use standard quantum control optimization techniques to 
find the T-long pulses which guarantee higher accuracy. 
The connection between OC and the QSL emerges at the 
time duration T < Tqsl for which OC fails to converge. 

Specifically, given an input state |V'(0)) and a Hamilto- 
nian H{t) that depends on the set of time-dependent con- 
trol functions x{t) = {xi{t), X2{t), ■ ■ ■ ,Xk{t)}, we shall 
employ the Krotov algorithm 0] to determine the optimal 
Xopt{t) that minimizes the infidelity X = l — \{tp{T)\ipG)\'^, 
which measures the distance between the target state 
IV'g) and the time-T evolved state \ip{T)) under H{t). 
The Xopt{t) are constructed iteratively, starting from 
some initial guess functions Xga{t). We then analyze the 
performance of the process as a function of T and show 
that the method is able to produce infidelities arbitrarily 
close to zero only above a certain threshold Tqsl, which 
we compare with the dynamical bounds that affect the 
system. We found a good agreement between these (in 
principle) independent quantities, meaning that the ef- 
fectiveness of our control pulses is only limited by the 
dynamical bounds of the system. Considering the lim- 
ited set of controls we allow in the problem, and the fact 
that our initial equations are not meant to optimize T , 
this is a rather remarkable fact that suggests that OC is 
a possible candidate for an operational characterization 
of the QSLs of complex systems. 

Even though our findings have been obtained in several 
different contexts, including for instance ordered Ising 
and Lipkin-Meshkov-Glick models, for the sake of clarity 
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FIG. 1: (Color online) Infidelity X versus number of itera- 
tions n of the Krotov algorithm 0] for different values of T 
(in units of /i/a;) for r(T)/aJ = 500. The dashed line is the 
estimated QSL (Tqsl = 1.5688) while the dot-dashed line is 
an exponential fit. Inset: Second derivative of the infidelity 
logarithm with respect to the logarithm of n. 

and for demonstration of the generality of the argument, 
in this Letter we shall focus on two paradigmatic exam- 
ples: the Landau-Zener (LZ) model [ll|, and the trans- 
fer of information along a chain of coupled spins with 
Hcisenberg interactions. The former case constitutes a 
basic step for the control of complex many-body systems, 
whose evolution, for finite size systems, is in many cases a 
cascade of LZ transitions [l^ . Adiabatic quantum com- 
putation is known to be limited by avoided crossings 
in the time-dependent system Hamiltonian and by our in- 
ability to avoid excitation of the system. The spin-chain 
case is instead related to one of the central requirements 
for the construction of circuit-model quantum comput- 
ers: an infrastructure that can rapidly and accurately 
transport qubit states between sites |14| . 

Landau-Zener model - The first example we consider is 
the paradigmatic case of the passage through an avoided 
level crossing 

W)] _r(,)), (1) 

in which V{t) is the control function that we shall op- 
timize through the Krotov algorithm. We start the 
evolution by preparing the system in the instantaneous 
ground state of -ff[r(0)] and we assume as our tar- 
get the ground state of H[T{T)], with r(r) = -r(0) 
(lim|r(o)|-,oo(^(O)|'0G) = 0). As an initial guess Tgsit) 
for the control we follow Ref. [l^. Here on the basis 
of the adiabatic theorem [l^ the control pulse T{t) was 
selected through a diflFerential equation F = 7 G^(F), 
where G(F) = -\- F^ is the instantaneous energy 

gap of the Hamiltonian ([T]), while 7 = [arctan(F(T)/a;) — 
arctan(F(0)/w)]/4Ta;. Starting from the Tgs{t) defined 
above, we run the OC algorithm for various values of the 
total time T. The results are reported in Fig. [Hby plot- 
ting the infidelity T as a function of the iterations n of 
the algorithm. When T < Tqsl ~ 1.5688, the infidelity X 
does not converge to zero, being its curvature asymptoti- 



FIG. 2: (Color online) Comparison between the time indepen- 
dent estimate (dashed line) and the second derivative criterion 
(circles, for -r(0)/cj = 5,10,20,100,500,5000 from right to 
left) for Tqsl for the LZ model. 

cally fiat. On the contrary, by progressively increasing T 
towards and above Tqsl, the curvature changes sign and 
the infidelity in the large iteration limit decreases expo- 
nentially, as confirmed by the fit in Fig. [TJ In the inset 
of Fig. [TJ data for the second derivative of the infidelity 
logarithm with respect to the logarithm of n for differ- 
ent T are shown: the derivative starts to cross the zero 
line for T w 1.58, and for T > Tqsl it clearly becomes 
negative. These findings are reflected by the study of the 
pulse shape of the optimization process (data not shown) . 
For T < Tqsl, the pulse develops a peak which grows in- 
definitely by increasing n and the control seems unable 
to converge towards an optimal shape. On the contrary, 
when T > Tqsl, after a certain number of iterations, the 
shape becomes stable, and only small corrections of the 
order of the infidelity take place. Remarkably, the pe- 
culiar feature of the initial guess Fgs (t) of being almost 
constantly zero for most of the central part of the evolu- 
tion is preserved by the recursive optimization of OC 0], 
suggesting that for this simple model an estimation of a 
finite resource QSL bound Tqsl for T can be deduced by 
a time- independent formula, assuming Hq = H[T = 0] as 
the Hamiltonian. In other words, for most of the evolu- 
tion time the dynamics can be effectively described by a 
time-independent Hamiltonian, which we can use to an- 
alytically estimate the QSL. This can be quantified with 
the Bhattacharyya bound Q , yielding 

Tqsl ^ AEq^ arccos |(V'(0)|Vg>| , (2) 

where AEq is the energy variance of Hq calculated on the 
initial state |i/'(0)), i.e. AT^o = [w^ - Auj'^/G^{T{0))]^/^. 
This approach has the advantage of providing a bound 
for T that is independent of the effective shape of the 
selected pulse. Finally, in Fig. [2] we show a comparison 
between the estimate Tqsl through the second derivative 
of the infidelity and the theoretical time-independent es- 
timate ^ for various F(0)/w ratios. We stress that here 
we have no fitting parameters. The excellent agreement 
shows that the OC efficiency is ultimately set by the dy- 
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namical bound of Eq. ([2]) . 

Quantum state transfer ■ 



We now apply our analysis 
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to a scheme for information transfer in a spin chain 
In this context a notion of QSL for spin chains was in- 
troduced in Ref.[l3l where the velocity of the informa- 
tion propagation was optimized with respect to the con- 
stant interactions among spins in the chain. This is of 
course quite different from the approach we introduce 
here where the couplings are given and the information 
transfer is sped up by using properly tailored external 
pulses. The model consists of a one-dimensional Hciseu- 
berg spin chain of length TV described by the Hamiltonian 

^ ^ , v^AT-l C{t) ( ,/,^^2 , 
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where a = {<Jx,'^y,'^z) are the Pauli spin matrices. J is 
the coupling strength between nearest-neighbour spins, 
C{t) is the relative strength of an external parabolic 
magnetic potential, and d{t) represents the position of 
the external potential minimum at time t (in units of 
h/J) along the axis of propagation Due to the 

conservation of the z component of the magnetization, 
we can restrict our analysis to the sector with a sin- 
gle spin-up only, so that a general state is described by 
1^(0) — X]m=o '^ml'^)' where |m) represents the state 
where the m-th site has its spin pointing up, and all other 
sites have spins pointing down. The states {\m)}^Z^ 
form a complete orthonormal basis for our Hilbert sub- 
space. Our goal is to evolve the initial state |V'(0)) = |0) 
to the final state \iI)g) = 1^ ^ 1): i-c to transport a spin 
up state from the first site to the final site of the chain. In 
Ref. [ill this was achieved by invoking the adiabatic ap- 
proximation, which relies on the fact that slowly moving 
the parabolic magnetic potential along the chain allows 
the spin-up to migrate from the leftmost site to adjacent 
sites via a nearest-neighbour swapping, while interac- 
tions between sites far from the field minimum are frozen. 
Specifically, the transfer was obtained by assuming pulses 
of the form C{t) = Co and d{t) = t{N - l)/r, and by 
working in a regime of large (ideally infinite) transmission 
time T. In our approach, we shall use instead the Kro- 
tov method to find the OC functions C{t) and d{t) using 
pulses similar to those of Ref. as an initial guess: this 
allows us to shorten the transfer time beyond what is al- 
lowed by the adiabatic regime. Once again, we optimize 
the controls for different values of T in an effort to iden- 
tify the minimum transfer time allowed by the selected 
controlling Hamiltonian. For each selected T the opti- 
mization algorithm stops either after a certain number n 
of iterations (of the order of 10^) or when the infidelity I 
reaches a certain fixed target threshold I* . The results 
obtained are shown in Fig. [3] and resemble those we have 
seen in the LZ model. In particular, as shown in Fig. [3] 
(left) the infidelity appears to converge to zero only for 
values of T that are above a certain critical time Tqsl- 
On the contrary, for T < Tqsl the convergence of the 
infidelity slows down, providing numerical evidence of a 
non-zero asymptote for n — > oo. In contrast with the 




FIG. 3: Left: Infidelity as a function of the number of it- 
erations n for a spin chain of length A'' — 101 and different 
durations T (in units of h/J). Right: Minimum infidelity T 
for chain length A'' and the (rescaled) transfer time (T — b) /N 
{b ~ 1.53 is the y-axis intercept of the function rQSL(A'^) re- 
ported in the inset), showing the expected linear behaviour 
of Tqsl with the size A'^. The green line follows from the es- 
timate Tqsl of the QSL time obtained by choosing the time 
T{N) at which the infidelity reaches the value T* = 5 x 10~^ 
for n = 10^ 



LZ case however, the dependence of 2 with the iteration 
number n is now less regular, refiecting the fact that the 
spin-chain dynamics is more complex than for the LZ 
model. Consequently, for the present model the sign of 
the second derivative cannot be used as a reliable signa- 
ture of Tqsl. Nonetheless, a numerical estimate Tqg^ 
for such quantity can been obtained by considering the 
smallest time T which allows us to achieve the target 
infidelity threshold X* in a fixed number of algorithm it- 
erations n (the result does not depend significantly on 
the value of I* and n). Apart from the case of small 
N, where boundary effects are more pronounced, the re- 
sulting TqgL appears to have a linear dependence on the 
chain length N - see Fig. [3] (right). 

For a comparison with an independent theoretical 
estimate of Tqsl, we cannot directly use the Bhat- 
tacharyya bound since in this case we are not al- 
lowed to treat the Hamiltonian as approximately time- 
independent. Nevertheless, a bound on the minimal 
transferring time can be obtained by considering the 
mean energy spread, obtained by averaging the instanta- 
neous energy spread of the system of the time-dependent 
Hamiltonian H{t) over the time evolution [0,T]. We de- 
fine this by A£a = 5^ AExit) dt, where AEx(t) = 
\J ((/)|[i/(t) — E\{t)Y\(^) is the energy spread and E\ (t) = 
{(f)\H {t)\(t)) is computed on the state \(f>), and A = 1,2 
labels different choices of For A = 1, we fol- 

low Ref. and take to be either the initial state 
\ip{0)) or the target state H'g), whichever results in the 
smaller A£i. For A = 2, we choose \(j>) = \ip{t)), which 
means we effectively divide the total evolution time T 
into smaller intervals dt over which H{t) can be assumed 
to be constant, and apply to each of them the Bhat- 
tacharyya bound 0] for timc-indcpendcnt Hamiltonians. 
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FIG. 4: Left: AE2{t) for a spin-chain of length TV = 101, 
and different total times T — 68.97 (solid green) and T = 
56.66 ~ Tqsl (dashed blue) (the times have been rescaled to 
make them fit on the same axis). The central portion of the 
plot has been enlarged in the inset. Right: '^2A£"i^^ (ors^nge), 

^2As^'' (I'sd circles), and Tqsl (green squares) versus chain 
lengtli N. 



obtained (for Ising coupling) in Rcf. [l^. Interestingly 
enough, however, our method only uses single-site local 
pulses while in Ref. [l^l this was achieved using global 
pulses that operate jointly on the whole chain. 

In summary, we have demonstrated that there are fun- 
damental constraints governing the efficiency of Krotov 
quantum OC algorithm dictated by the maximum speed 
at which a quantum state can evolve in time. These re- 
sults provide a further link between control theory and 
quantum dynamics. 

We thank A. Carlini for discussions, and the bw- 
GRiD (http://www .bw-grid.de) for the computational 
resources. We acknowledge financial support by 
SFB/TRR21, and the EU under the contracts MRTN- 
CT-2006-035369 (EMALI), IP-EUROSQIP, and IP- 
SCALA. 



The bound on the minimum transfer time is then given 
by rQSL/(A^— 1) > max{ ^^^^ ; }, which needs to be 
satisfied by any H{t) that brings |V'(0)) to the target state 
lipc) in ^ time T. Since we have taken the time average of 
the bound, we interpret the quantum speed limit as de- 
scribing the minimum transfer time 'pcr-site'. In Fig. [J] 
(left), we report Ai?2(t) as a function of time for two dif- 
ferent total transfer times. The energy spread is almost 
constant save near the final time, where large oscillations 
are present, corresponding to the deceleration of the spin- 
wave. The picture that arises is that OC finds a solution 
Xoptit) that initially accelerates the spin excitation, and 
then transfers it with constant velocity up to the end of 
the chain, where deceleration occurs. The average en- 
ergy spreading AE2 (t) is larger for smaller time transfer 
T allowing for a higher average excitation velocity. We 
finally compare the estimated optimal time Igg^ from 
Fig- [3] (right) with the analytical estimate of lihe QSL 
given above. The results are reported in FigH] (right) 
where the two quantities are compared: both estimators 
are linearly dependent on N, but the numerical results 
show an improvement over the theoretical prediction by 
a factor of 77 « 3, which we attribute to the difficulty in 
formulating the quantum speed limit for our many-body 
problem. One can think of the optimal transfer of the 
excitation as being facilitated by a cascade of effective 
swaps. As before, the agreement of scaling of the two 
results shows that even in the presence of additional con- 
straints the OC reaches the ultimate dynamical bounds 
set by quantum mechanics. Indeed, we achieved an im- 
provement of the transfer time with respect to Ref. [l8| 
of up to two orders of magnitude. It is also worth men- 
tioning that we achieved a transfer time faster than that 
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